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[i] To investigate ion density depletion along magnetic field lines, we compare in 
situ-measured ion density fluctuations as seen from C/NOFS and compare them to the 
field-line-integrated depletion of the whole bubble as inferred from electric field 
measurements. Results show that, within C/NOFS’ range, local measurement of the 
normalized density depletion, A/i/zio, near the apex may be far less than at other points on 
the same field line. We argue that the distribution of AnJn 0 is a weighted distribution 
concentrated at latitudes of the Appleton anomalies and becomes more heavily weighted 
the closer the field-aligned bubble rises to the peak of the anomalies. A three-dimensional 
simulation of an ionospheric bubble verifies our arguments. 
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I. Introduction 

[ 2 ] Since sprcad-F irregularities are field-aligned, the 
spread-F problem is often simplified by integrating three- 
dimensional (3D) physics along magnetic field lines into a 
two-dimensional problem, as first prescribed by Perkins 
[1973] for midlutiludcs and by G. Hacrendel (Theory of 
equatorial spread F, preprint, Max-Pianck-lnstitut fur Extra- 
icncstrischc Physik, Munich. Germany, 1973). However, to 
locate where scintillation will lutiludinally occur, the distri- 
bution of irregularity along a field line is important. The 
problem is obfuscated by a lack of spread-F measurements 
along field lines. No radar is capable of the required latitu- 
dinal range. Satellite orbits cannot follow along field lines. 
Hence, we arc limited to local measurements of spread F and 
lo only a one-dimensional cut across a bubble. However, 
because of the high conductivity parallel to Earth's magnetic 
field. B. field lines are very accurately cquipotential for 
scales greater than 1 km [Kelley, 2009], which is why spread 
F is field-aligned in the first place. From a single local electric 
field measurement, we arc able to infer the electric field at 
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any poinl along the field line, including those aligned with 
spread F. The electric field mapping along magnetic dipole 
field lines is derived in Appendix A. We later show that, from 
a single local electric field measurement, we can also infer 
the field-linc-intcgrated ion density depletion. We analyze 
and interpret in situ measurements of ion density irregulari- 
ties and compare them to infrared ficld-line-integratcd 
depictions from electric field measurements as observed by 
the Communications/Navigation Outage Forecasting System 
(C/NOFS) satellite. 

[j] Developed by the Air Force Research Laboratory. 
C/NOFS was launched to investigate the F region of the 
equatorial ionosphere. The vehicle orbits at an inclination of 
13 degrees in reference to the geographic equatorial plane. 
The satellite orbit has a perigee of 400 km and an apogee of 
850 km. Part of the satellite's payload is a Planar Langmuir 
Probe (PLP) [Roddy el al., 2010], which measures ion 
density. C/NOFS also measures electric fields with the 
Vector Electric Field Instrument (VEFI) [Pfaff el al., 2010], 
With both PLP and VEFI. there is an abundance of simulta- 
neous in situ ion density and electric field measurements. The 
data presented in this paper were measured from May 2008 
through October 2009. Using PLP to measure the local ion 
density depletion at a local point and VEFI to infer the 
depletion along the whole field line, we explore the distri- 
bution of density depiction along field lines. To infer the 
depiction along an entire field line, we must first understand 
the dynamics of a bubble in three dimensions. 

2. The Ion Density Perturbation Growth 
of Irregularities Along Field Lines 

[4] Throughout this paper, we work in the frame of view 
moving with the E x B drift of the background electric field. 
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E 0 . As derived in Appendix B. the current equation in this 
frame is 


J = <r -AE + 


nils' x B Vp x B 


<D 


where J is the current density, a is the conductivity tensor, 
E is the electric field, AE is a perturbation electric field 
(AE = E - E 0 ), n is ion density, M is ion mass. B is the 
magnetic field, and p is pressure. Note in this form that the 
background electric field is found in the effective gravity jj' 5 
jj — I'faEo X B JB 2 where g is gravity, r-m is the ion-neutral 
collision frequency, and Eg is the background electric field. 
Assuming that Vp is vertical and that gravity is anti- vertical 
and they only drive zonal currents and zonal polarization 
electric fields, we restrict our focus to the zonal components 
of current equation (1): 


7,=o,A£,+^cosD-^|cosD, 


< 2 > 


where On is the Pedersen conductivity, A E, is the zonal 
electric field perturbation, and D is the dip angle (the angle 
between the magnetic field and the horizon). Throughout 
the paper, we refer to .v as being in the magnetic zonal 
direction, y as being parallel to the magnetic field, and z as 
being perpendicular to the magnetic field and in the meridional 
plane. The problem with (2) alone is that we must recognize 
that all of the points along a field line are approximately 
equipotential, as discussed earlier, and arc thus effectively 
tied together. Every point along a field line has the same net 
zonal current as any other point on the line. Note that current 
observed at any specific point is not necessarily produced at 
that point locally. For instance, electric field-induced currents 
will flow through the paths of least resistance and not nec- 
essarily where the electric field is strongest. Regardless of 
how electric fields map along magnetic field lines, most of 
the current will flow through where the Pedersen conductivity 
is greatest. A more comprehensive current equation must be 
formulated. 

[s] Take the circuit analogy of several current sources tied 
together in parallel. The output current of the parallel sources 
is the sum of all the sources. In our case, the current sources 
arc driven by electric fields, gravity, and pressure gradients 
tied together in parallel by magnetic field lines. The total 
zonal current driven by any of the three drivers is a field-line 
integral of its current contribution. Dividing this integral by 
the length of the line yields the effective zonal current den- 
sity. Using the electric field mapping for a dipole magnetic 
field (Al), the effective zonal current driven by electric fields 
is formulated as 


h 


<3> 


where AE.)*-,, is the electric field perturbation at the mag- 
netic equator, 0 is the magnetic latitude, and the integration in 
the y direction is along the magnetic field line. Turning to the 
other current contributions, we note that Perkins [ 1 973] treats 
both gravity and prcssure-gradient-driven currents as a 
varying value along the field line. However, similar to the 


formulation of Pedersen currents, both gravity -driven cur- 
rents and pressure-driven currents arc also tied together by 
magnetic field lines. The effective zonal currents should 
not differ from point to point along equipotential field 
lines. All points on a field line must be subjected to the same 
net zonal current, regardless of how much current is pro- 
duced locally. Similar to the formulation of (3), we formu- 
late the effective current density driven by gravity and 
pressure gradients in equations (4) and (5), respectively, as 
field-line integrals divided by the length of the field line: 
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where D is the dip angle, the angle between the magnetic 
field and the horizon. Combining all of the effective current 
densities yields 


r,= 
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> 4 - / ^ cosDdv 

-AE,\„+ J B 


7 * 


/ ** 


/is 00 ^ 

“7^ 


(6) 


where the first, second, and third terms are the effective 
current densities driven by electric field, gravity, and pres- 
sure gradients, respectively. The new current equation (6) is 
similar to the equation of the locally produced current (2) 
except that all of the points along a field line are properly 
tied together and are actually common values to all points on 
equipotential field lines. 

[6] To satisfy the quasi-neutrality of plasma, the current 
divergence must be zero ( V • J = 0). Assume the plane wave 
solutions AE, = 6E^ U " * >l and A n = in e* 1 "' ta> where k is 
the zonal wave number and An is the ion density perturba- 
tion from the ambient density, n 0 . The corresponding pres- 
sure is p - p 0 +■ k B (T, + T,)An where p 0 is the ambient 
pressure, k B is the Boltzmann constant, 7, is the ion tem- 
perature, and T, is the electron temperature. Taking the 
divergence of (6), setting it to zero, using the plane wave 
solutions, dropping second-order terms, and using the electric 
field mapping for dipole magnetic fields (Al), the lineariza- 
tion of the zonal current divergence equation yields 


*£,| # -o = - 


/ 


Mg’ cos D 


indy 


/ -^h dy 


AE, = — 


./ 


Mg' cos Z> 
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Andy 


(7) 


where 6E , |^-o is 6E„ evaluated at the magnetic equator, 
and we define an effective Pedersen conductivity as Op — 
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cos 5 ( 0) • / {dy ■ oj cos 5 Oyjdy. We neglect the second-order 
term involving the change in a p due to density perturbations 
for two reasons. For cases where An and A E, arc small 
compared to the ambient density and ambient zonal electric 
field, respectively, the second-order term is negligible. Sec- 
ond-order terms arc often neglected in linear analyses. Even 
if An and A E, are not small, the change in the effective 
Pedersen conductivity, is negligible. Depletions arc 
found in the F region and, while they affect o p in the F 
region, they will not significantly affect a p at lower altitudes. 
Since the effective Pedersen conductivity is dominated by 
a. in the minimally affected lower altitudes, neglecting the 
changes of a p due to density perturbations is justified. 
Equation (7) describes the polarization electric field that 
develops to relieve any current divergence due to density 
perturbations along a field line. Note that since d(A n)'dz is 
zero and the zonal divergence of cptf'dz is zero, pressure 
gradients have no zonal current divergence in our analysis 
and play no role in (7). A circuit analogy to AE, as for- 
mulated by (7) is the voltage that results from several current 
sources and conductors in parallel with one another. TTic 
voltage across such a circuit w ould be the total current sources 
divided by the total conductance. In our case, the divergences 
of gravity-driven currents are current sources and the Pedersen 
conductivities are conductors, both of which arc tied together 
in parallel by magnetic field lines. 

[7] To formulate the ion density, perturbation growth of a 
bubble at a given point, we turn to the continuity equation in 
a frame of view moving with the gravitational and diamag- 
netic drift in addition to the frame moving with the E X B 
drift of the background electric field. Bubbles move with 
both gravitational and diamagnetic drifts and remain zonally 
fixed in this frame of view. The remaining velocity to con- 
sider is the E X B drift of AE,. Assuming the plasma is 
incompressible (V • v = 0) and neglecting recombination, 
which we assume is much slower than the processes of 
spread F, the continuity equation becomes 


J + V(nv) = 0, 

dn _ cm AE, 
di~ 


(8) 


which describes the growth rate of the density perturbation 
at any point in a bubble. From (7), AE, is the opposite sign 
of An. Hence, density perturbations grow’ where cn^dz is 
positive. In other words, density perturbations grow below the 
F peak. For a sense of the variation of the density perturbation 
growth rate (8) along a field line, we model the varying B and 
AE, along a field line. Modeling B as a magnetic dipole such 

that B= (\/l + 3 sin 2 0/cos 6 0) ■ where is the 

magnetic field at the equator and using the zonal electric 
field-mapping equation (Al), we can recast (8) as 


dn 

di 


cos ! 0 dn,, AE, 

\/l + 3 sin 2 0 s l.-. 


fdno AE,\„ 
v/l+3sinW * «L, 


(9) 


[s] Note that A£,|^ and arc global values shared by 
all points on the same field line whereas 0 and cn<ydz are 
local values. How’ much the density perturbation at a point 
along a field line grows relative to other points on the field 
line is determined by its magnetic latitude and the local 
gradient of the background ion density. It is apparent from 
the cos J 0/\/l + 3 sin 2 0 factor that the greater the distance 
from the equator, the greater the local density perturbation 
growth attenuation. The attenuation of irregularities away 
from the equator can be seen in the work by Dao el al. 
[2011, Figure 2] where plasma irregularities are shown to 
straddle the magnetic equator, as observed by C/NOFS. 

3. Inferring Field-Line-Integrated Ion Density 
Depletion From Electric Fields 

[9] From equations (7)-f9) describing the dynamics of 
bubble growth, we can infer the ficld-line-integratcd deple- 
tion from a single measurement of AE, at any point on the 
field line. In equation (7) describing the polarization electric 
field induced by a bubble, we believe the integral in the 
numerator is dominated at the F peak. As suggested by the 
growth rates in equations (8) and (9), An maximizes near 
the F peak. Also, considering that most of the gravity- 
driven currents exist at the F peak where the ion density is 
greatest most of the current divergence must happen near 
the F peak. We thus approximate equation (7) as 


( 10 ) 


where we define 


[ Mg' cosD 

Mg' cosD\ _ f b dy 
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01 ) 
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dy 


as a weighted average, concentrated where we assume the 
majority of An is located: near the F peak. The variable q is a 
tunable parameter controlling how much weight to empha- 
size at the F peak. Modeling B and the field lines with the 
International Geomagnetic Reference Field (1GRF), M from 
the International Reference Ionosphere (IRI) model [Bilitza 
and Reinisch, 2008], and a p and from both the IRI and 
the mass spectrometer incoherent scatter model (MS1S), we 
can infer the total depletion of a bubble along a field line from 
VEFI after re-arranging equation (10) as 



[10] The field-line conductivity is still greater than the 
Pedersen conductivity into E region altitudes, even at night. 
Hence, when inferring AN, wx integrate as low as 90 km. 
We do not consider anything below 90 km altitude because 
the increased ion-neutral collisions reduce the Pedersen 
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Figure 1. Average Anln 0 :AN/N 0 from May 2008 to 
October 2009 where A N/N 0 is inferred from A E,. The 
dotted-dashed lines represent dipole magnetic field lines. 
The higher the apex height of each field line, the smaller 
&n/n 0 :&N/N 0 is, indicating a much greater A n/n 0 some- 
where farther down the intersecting field line. 


conductivity significantly and we expect minimal zonal 
cunent below this altitude. Note that because we do integrate 
to such low altitudes, one cannot approximate the ion gyro 
frequency to be much greater than the ion-neutral collision 
when computing the local Pedersen conductivity, which is 
Op = futfilhUyh, + f*i) where H, is the ion gyro frequency. 

4. Data Representation: C/NOFS Ion Density 
and Electric Field Measurements 

[n] From the C/NOFS’ PLP, we derive the ambient den- 
sity. no- as 60-s (~4S0 km spatially) averages of measured 
density. Subtracting n 0 from the density measurements yields 
the density deviation. An. We normalize the deviation by the 
background density as An/no. We determine A£, from VEFI 
similarly as An, deriving the ambient zonal clcctnc field as 
60-s averages of the zonal electric field measurements and 
then subtracting it from the same measurements. Using 
equation (12), we then infer AY from A E,. We also compute 
No • fnody from 1R1 and normalize the ficld-linc-intcgratcd 
depiction as A N/N&. We use the normalized parameter 
AN/No as opposed to the absolute AY because if the den- 
sity profiles from 1R1 are off by a factor when computing 
the numerator, this (actor is roughly canceled out by the 
denominator. We now focus our attention on the ratio 
An/n 0 :AVW 0 . which relates the local A n/n 0 to that of the 
whole field line. The ratio indicates to what extent the local 
An/no exceeds that of the whole field line (An>io:AA'/.V) > 1) 
or is below that of the whole bubble (An/n 0 :AAW 0 < 1). 
Note that An/n 0 :AN/N 0 is highly sensitise to noise as AY/Y 0 
approaches zero. To remedy this issue in our statistical study, 
we compute the 60-s mean of the magnitude of A/i/ho and 
AY/Y 0 separately and present An/n 0 -.AN/N 0 as the ratio 
between these two mean magnitudes. We restrict our focus to 
irregularities, so in our data presentation, we only use data 
where either of the mean magnitudes of An/n 0 or AV.Y 0 


exceeded 10%. TTie average A/i//i 0 : AYW 0 is binned by 
magnetic latitude and altitude, as shown in Figure 1. Die 
dotted-dashed lines trace dipole magnetic Geld lines. 'Die 
ratio An'n 0 :AN/N 0 is nearly symmetric across the magnetic 
equator. We attribute the slight asymmetry to the asymmetry 
of Earth's magnetic field. 

5. Discussion 

[ 12 ] As evident from equations (7) and (9), determining 
the distribution of ion density along the field line can be 
complex, beginning from an initial perturbation and then 
depending on the time evolution of (7) and (9). The com- 
plexity results from the local growth rate of An being 
dependent on the collective An along the whole field line, as 
dictated by (7), while the growth rate of An all along the 
field line also depends on unique local conditions. Further 
complicating the problem, the rise velocity of the bubble 
varies as a function of altitude, bubble size, and the time- 
varying background electric field [McDaniel and liysell, 
1997]. These considerations aside, we consider another 
important factor in the distribution of ion density irregularity 
along a field line: the latitudinal variation of plasma density 
profiles. To illustrate, we make a zero-order approximation 
that the ratio of AEJ^/B,^ to the rise velocity ( V = dzJdt) 
is constant, the cos J 0/(1 +3 sin 2 0) factor is approximately 
unity, and the gradient of n 0 is roughly in the z direction. 
Equation (9) can then be approximated as 

<,3 > 



-20 -10 0 10 20 


magnetic taUude (flog) 

Figure 2. Schematic of a bubble rising through the Apple- 
ton anomalies where gray shading indicates the ion density 
of the anomalies. Die bubble locally grows most rapidly 
where the gradients in the z direction arc the greatest beneath 
the anomalies' peaks. The resulting An/n 0 distribution con- 
centrates near the anomalies. 
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Figure 3. Snapshot of 3D bubble simulation, (top) An equatorial slice of the bubble's ion density in units 
of log(cm *). The dotted-dashed lines indicate the meridional slices that arc examined in their respective 
columns, (middle) The ion density of the mentioned meridional slices, (bottom) The computed An/ngi 
AN /No of the bubble. The part of the bubble rising closer to the Appleton anomalies yields a heavier 
distribution of An/n 0 near the anomalies. 


[u] Equation (13) suggests that An is roughly propor- 
tional to n 0 . If the density profile was latitudinally indepen- 
dent along the magnetic field line, the distribution of 
normalized density depiction. An/no, will be approximately 
uniform after dividing each side of (13) by n 0 . This would be 
a tremendous coincidence, allowing us to infer An at any 
point along the same field line from a single measurement of 
An/n 0 by simply inferring n 0 at the point of interest with an 
ion density model. However, density profiles do vary with 
latitude, especially in the presence of the Appleton anoma- 
lies. Consider the schematic of a bubble rising through the 
anomalies shown in Figure 2 where the gray shading repre- 
sents the ion density of the anomalies. Before approaching 
anomalies where the density profiles vary less significantly 


with latitude, we can expect a more uniform distribution of 
An/n 0 throughout the bubble. However, as a bubble rises, 
the local points approaching the anomalies experience greater 
density /-gradients than other points along the bubble. Con- 
sequently, the density perturbation grows rapidly at the 
points near the anomalies and the resulting An!n 0 distribution 
becomes heavily weighted there as well. Because most of the 
depiction occurs near the anomalies. An /n& AN/N 0 increases 
near the anomalies while decreasing at other points such as 
the apex (where the ONOFS observations presented in 
Figure 1 were observed). 

[m] During 2008 and 2009, when the data presented were 
collected, the solar flux was at a minimum at uround 60 solar 
flux units. Under such conditions, the F peak and the 
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Applclon anomalies were aboul 300-350 km high (below 
C/NOFS) and the apex altitude (the height of the apex of 
the intersecting field line) of the Appleton anomalies was 
about 800-900 km high (the upper limit of C'NOFS' range). 
Within C/NOFS’ range (400-850 km), the observed irregu- 
larities with higher apex altitudes were ones that have risen 
closer to the Appleton anomalies. We expect the irregularities 
with higher apex altitudes to have a Ant'n 0 distribution 
increasingly concentrated near the anomalies, resulting in 
decreasing An/n 0 :AN/N 0 at CNOFS’ point of observation. 
As seen in Figure I, An/n^AN/No actually decreases with 
increasing apex altitudes. From a more careful inspection of 
An/n 0 :AN/N 0 along any given field line, An/n 0 :AN/N 0 
increases farther away from the equator. We attribute this 
phenomenon to increasing density /-gradients away from the 
equator, approaching the latitudes of the anomalies. Unfor- 
tunately, C/NOFS* range does not extend to the anomalies at 
any time where we expect An/n 0 :AN/N 0 > I. We turn to 
numerical 3D modeling to check the consistency of our 
arguments. 

[is] We present the use of a 3D ionospheric model out- 
lined by Retterer [2005, 2010] to calculate both global 
ambient ionospheric structure and 3D bubble simulation. 
The ambient model solves the continuity equation for 
plasma density in terms of the processes of production, loss, 
and transport of plasma density along and perpendicular to 
geomagnetic flux tubes in the vicinity of the geomagnetic 
equator. It requires specification of the solar UV flux, neutral 
densities, winds, and temperatures within the thermosphere, 
along with plasma drift velocity and temperature. For this 
study, empirical models of these parameters were used: the 
Mass Spectrometer and Incoherent Scatter (MSIS) model 
[I/edin, 1987] for neutral densities and temperatures, the 
Horizontal Wind Model (HWM) [Hedin el al., 1991] for 
neutral winds, the IlmtereRger el al. [1981] model for the 
solar ultraviolet spectrum, the Scherliess and Fejer [1999] 
model for plasma drifts, and the Gulyaeva and Tilheridge 
[2006] model for the plasma temperature. We specify the 
F10.7 to be 60 solar flux units in the simulation presented in 
this paper. After the ambient conditions are computed, a 3D 
bubble simulation is done as described by Retterer [2010]. 

[is] A sinusoidal perturbation with a wavelength of 1200 
km was initiated and we took a snapshot of the bubble as it 
rose to an altitude of 600 km, as shown in Figure 3. Figure 3 
(top) is an equatorial slice of the bubble’s ion density in units 
of log(cm 3 ). The dotted-dashed lines indicate the meridio- 
nal slices that arc examined in their respective columns. 
Figure 3 (middle) presents the ion density of the mentioned 
meridional slices. We compute Anlno'.AN/No of the bubble 
in Figure 3 (bottom). In Figure 3 (left), a meridional cut of 
the bubble far below the anomalies, the An/n 0 distribution is 
nearly uniform with a An/n 0 :AN/N 0 near unity, indicating 
that An is roughly proportional to n 0 . However, in Figure 3 
(right), we show a slice of the bubble rising closer to the 
anomalies. As expected, the w eight of the A n/n 0 distribution 
lies heavier where the bubble intersects with the Appleton 
anomalies. Consequentially, An/n 0 :AN/N 0 is significantly 
greater than one at the Appleton anomalies and significantly 
less than one about the equator where the bubble docs not 
pass though the anomalies. This is consistent with C/NOFS 
data representation and supports our argument of a heavy 


distribution of Anin 0 near the Appleton anomalies as bub- 
bles approach them. 

6. Conclusion 

[ 17 ] The distribution of an ion density bubble along a field 
can be quite complicated since it depends on the time history 
of the bubble rise, starting from the initial perturbation and 
continuing with varying local conditions that impose vary- 
ing local growth rates of An along the bubble. In the pres- 
ence of Appleton anomalies, Antn 0 cannot be approximated 
as being constant along field lines. The data presented in this 
paper suggest that any local measurement of An/n 0 near the 
apex of a field line is only a very small fraction of the total 
field-linc-intcgruted AN/N 0 , more so when the bubble 
intersects the Appleton anomalies. We note that a lot of 
variation existed in the data. This comes as no surprise, since 
longitudinal variability (realistic magnetic field. Appleton 
anomalies, tidal effects, etc.) as well as day-to-day variability, 
such as varying levels of solar flux, occurred during the span 
of 2008 to 2009. 

[is] An implication of this study is that local measure- 
ments of An/no or An at the equator cannot be trivially 
extended to higher latitudes in an effort to calculate scintil- 
lations at these latitudes. However, with equation (12), we 
show a means to estimate the density depletion of whole 
field lines from local zonal electric field perturbations. A E x . 
Equation (13) suggests that the distribution of An is similar 
to the distribution of no along the field line, so one may infer 
A n/n 0 or An at any point along a field line from a single 
measurement of AE, in conjunction with an ion density 
model. This calls for further investigation. We plan to test 
this approach using scintillation data from the —74° longi- 
tude chain [Valladares el al, 2004], Since C/NOFS’ launch, 
up until now, scintillations have been low at GPS frequen- 
cies, but as solar activity increases, this comparison may be 
possible. 

[ 19 ] Also, for increased solar activity, we hope to perform 
the same calculations as presented in this paper. The most 
significant change would be a higher F peak associated with 
higher solar flux. Under these conditions, C/NOFS may 
intersect the Appleton anomalies where we expect A n/n 0 : 
AA’/A'o to exceed unity. 

Appendix A: Electric Field Mapping 
for Magnetic Dipole 

[ 20 ] While magnetic field lines arc cquipotential, the dis- 
tance between field lines docs vary. Approximating the 
Earth’s magnetic field as a dipole magnetic field, r- Ro cos 2 0 
where r is the distance from the center of the Earth, Rq is 
the r of the field line at the magnetic equator, and 0 is the 
magnetic latitude. The zonal (x direction) electric field along 
a field line can be formulated as 


</«t> _ d 4> 

* ~~ ~dx ~ rcosOdX. 


E, = - 


E, = 


R 0 cos’&A. 


(Al) 


COS >0 
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where <I> is the potential. A is magnetic longitude, and £,1*-© 
is the zonal electric field at the equator. Using the magnetic 
dipole identity, sin D - dzl{r dO), the derivative of Ro with 
respect to magnetic latitude dRfJdO - 2 r sin fl/cos 3 0, and the 
dipole relationship sinD = 2sin OjJ 1 + 3 sin J 0 where D is 
the angle between the Earth's magnetic field and the hori- 
zon, the meridional (z direction) electric field along a field 
line is formulated as 


take x and z to be perpendicular to the magnetic field in the 
zonal and meridional directions, respectively. Subtracting 
the ion velocity from the electron velocity and multiplying the 
difference by the ion’clcctron density and electron charge, e, 
yields the current density equation in the frame of view 
moving with the E X B drift of the background electric field: 


J = a • AE + 


nMg x B 


Vp x B 


(B4) 


*/<D _ J«> 
d; rdO sin/) 


E, = J4>- 


E,= 


1 2/- sin# y'l 4- 3 gin 2 # 
r JRocos 3 # 2 sin# 

JO \/l +3 sin 2 # 

JR^ cos 3 # ’ 


E,= 


VI + 3 sin 2 # 

co?# 




(A2) 


where £j|*-c ' s *^ c meridional electric field at the equator. 
Both (Al) and (A2) are equivalent to expressions given by 
Mozer [19701, generalized for an arbitrary point on a field 
line instead of strictly where the field lines intersect the 
Earth's surface. 


where M is the ion mass and we assumed that k >• 1. The 
conductivity tensor, a, is 

/ Op 0 o k \ 

u = I 0 a 0 0 , (B5) 

\ -(7* 0 Op) 


where a p . <7j, a©, are the Pedersen, Hall, and parallel con- 
ductivities, respectively. 
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Appendix B: Current Density in a Frame Moving 
with E x B Drift of the Background Electric Field 

[ 21 ] The momentum equation in the neutrals frame of 
view is well known and can be written as 

itn^v = /w(E + v x B| + nmg - Vt> - nmv.v. (Bl) 

at 

where n is the ion or electron density, m is m ass, v is ion 
velocity, q is charge, E is the electric field, B is the magnetic 
field, g is gravity, V/> is the pressure gradient, and v„ is the 
collision frequency with neutrals. In the frame moving with 
the E x B drift of the background electric field, Eo. we make 
the transformation, v -* v + E„ X B ‘B 2 . Equation (Bl) 
becomes 


Eo x 


run — = nq( AE + vxB)4- nma - Vn - nmv, ( v 

at \ u' 

rtm — = nq( AE 4- v x B) + nmg - V/» — nmv.v, 




(B2) 


where we define the effective gravity, g' = g' + t>. E 0 x BIB 1 , 
as was done by Ott [1978]. Solving for v in steady state 
(d\!dt ~ 0) yields 



\ « J + 1 


0 



.( — AE 4- — - , (B3) 

,mv, v, wm> 


1 

«* + !/ 


where k is the gyro-to-collision frequency ratio (x - qBimv.) 
and we take y to be in the direction of the magnetic field. We 
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